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Recent advance in radioactive ion beam experiments has opened up a new field of research 
in nuclear physics (for a recent review, see, e.g., refs. [!ljj.j2|). These experiments have already 
provided useful information about the structure of unstable nuclei far from the stability 
valley. However, it has not been generally recognized that they are also useful for extracting 
information on the equation of state (EOS) of asymmetric nuclear matter, which has not been 
well determined and is important for understanding both the structure of unstable nuclei 
and the properties of neutron stars. In relativistic mean-field (RMF) theory, it has been 
recently shown that the nuclear symmetry energy affects significantly the binding energy 
and mis radii of neutron-rich nuclei |J. Also, the chemical composition, the evolution of 
lepton profiles, and the neutrino fluxes in neutron stars depend strongly on the nuclear 
symmetry energy [|]. More important but poorly known is the density dependence of the 
nuclear symmetry energy. Better knowledge on this quantity is required to understand both 
the matter radii of many neutron-rich isotopes, which has been found to increase faster than 
A 1//3 ||, and the central density of neutron-rich nuclei, which is lower than that of stable 
nuclei |J. In nuclear astrophysics, the density dependence of nuclear symmetry energy is 
crucial for understanding the supernova explosion scenarios and the cooling mechanisms of 
neutron stars [j]||. In particular, it determines the equilibrium concentration of protons in 
a neutron star. If the latter is larger than a critical value of about 15%, the direct URCA 
process can happen, and would then enhance neutrino emissions and the neutron star cooling 
rate @. 

Muller and Serot have recently shown that the asymmetric EOS has quite distinct new 
features compared to the symmetric one ||. In particular, in asymmetric nuclear mat¬ 
ter the liquid-gas phase transition is second order rather than first order as in symmetric 
nuclear matter. Furthermore, the instabilities that produce a liquid-gas phase separation 
in asymmetric nuclear matter arise from fluctuations in the neutron/proton concentration 
(chemical instability) instead from fluctuations in the baryon density (mechanical instabil¬ 
ity) as in symmetric nuclear matter. The latter may be related to the recently observed 
isospin-dependence of nuclear multifragmentation in heavy-ion collisions at intermediate 
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energies mm 

Presently available radioactive ion beam facilities and planned isospin laboratories pro¬ 
vide the unique opportunity to study experimentally the EOS of asymmetric nuclear matter. 
Tanihata has recently proposed to extract the EOS of asymmetric nuclear matter by studying 
the properties of neutron-rich nuclei, such as their density distributions, radii and the thick¬ 
ness of the neutron skin ||. In this Letter, we shall suggest a novel approach for studying 
the asymmetric part of the nuclear EOS, i.e., via the ratio of the number of preequilibrium 
neutrons to that of protons (i? n / p (fAi n ) = dN n /dN p ) from collisions of neutron-rich nuclei at 
intermediate energies. Using a transport model with explicit isospin degrees of freedom, we 
will show that the ratio R n / P (E kin ) is sensitive to the density dependence of the symmetry 
energy, but is almost independent of the compressibility of symmetric nuclear matter and 
the in-medium nucleon-nucleon cross sections. 

Theoretical studies (e.g., [fl],|T^-[T7|] ) have shown that the EOS of asymmetric nuclear 
matter can be approximately expressed as 

E(p,f3) = E(p,f3 = O) + S(p)0*, ( 1 ) 


where p = p n + p p is the baryon density; f3 = ( p n — p p )/{p p + p n ) is the relative neutron 
excess; and E(p,/3 = 0) is the energy per particle in symmetric nuclear matter. The bulk 
symmetry energy is devoted by S(p ) = E(p,/3 = 1) — E(p,/3 = 0). Its value S 0 = S(p 0 ) at 
normal nuclear matter density is known to be in the range of 27-36 MeV |L8|. In the non- 
relativistic Hartree-Fock theory (e.g., [|19|,^0[]) and the relativistic mean field (RMF) theory 
(e.g., |2lH27f ), the predicted values of S 0 are 27-38 MeV and 35-42 MeV, respectively. Also, 
the density dependence of the symmetry energy varies widely among theoretical studies. A 
p 1 / 3 dependence was obtained by Siemens using the Bethe-Goldstone theory for asymmetric 
nuclear matter |L2|, while the RMF theory predicts a linear dependence p8| , [29| . One of the 
most sophisticated calculations is that of Wiringa et al. using the variational many-body 
theory ||T5||. Different density dependences of S'(p) have been found depending on the nuclear 
forces used in the calculation. Typical results of these studies can be parameterized by |T6 
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(2) 


S(p) = (2 2/3 - 1)-E° f [u 2/ 3 - F(«)] + S 0 F(u), 


with F(u) having one of the following three forms 

pn 2u 2 
Fi(m) = 


1 + u 
F 2 (u) = u, 

F 3 (u ) = m 1 / 2 , 


(3) 


where u = p/po is the reduced baryon density and Ep is the Fermi energy. From eq. @ the 
contribution of nuclear interactions to the symmetry energy density can be obtained, i.e., 

w a (p,(3) = e a -pF(u){3 2 , (4) 

where e a = [{S 0 - (2 2 / 3 - 1)| E° F ] is the contribution of nuclear interactions to the bulk 
symmetry energy at normal nuclear matter density. The mean-field potentials for neutrons 
and protons due to the symmetry energy are then 

Ksy P) (/h P) = dw a (p, P)/dp n(p) . (5) 

To illustrate the magnitude of the symmetry potential, we show in Fig. srsyw) 
using the three forms of F(u) and So = 32 MeV. It is seen that the repulsive (attractive) 
mean-held for neutrons (protons) depends sensitively on the form of F(u), the neutron excess 
(3, and the baryon density p. In collisions of neutron-rich nuclei at intermediate energies, 
both (3 and p can be appreciable in a large space-time region where the isospin-dependent 
mean-fields, which are opposite in sign for neutrons and protons, are strong. This will affect 
differently the reaction dynamics of neutrons and protons, leading to possible differences in 
their yields and energy spectra. Since the magnitude of the asymmetric part of the nuclear 
EOS is small compared to the symmetric part in eq. (]!]), to extract S{p) requires observables 
which are sensitive to the asymmetric part but not the symmetric part of the nuclear EOS. 
Also, these observables should not depend strongly on other factors that affect the reaction 
dynamics, such as the in-medium nucleon-nucleon cross sections. In the following, we shall 
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demonstrate that the ratio R. n / P {E^ UI ) of preequilibrium neutrons to protons from collisions 
of neutron-rich nuclei meets these requirements. 

We shall use an isospin-dependent Boltzmann-Uehling-Uhlenbeck (BUU) transport 
model (e.g., |30|-|33|). The proton and neutron densities calculated from the nonlinear rel¬ 
ativistic mean-held (RMF) theory [2HI54] are used as inputs to initialize the BUU model 
35| , |36[| . The isospin dependence is included in the dynamics through nucleon-nucleon col¬ 


lisions by using isospin-dependent cross sections and Pauli blocking factors, the symmetry 
potential V£^\p,f3) and Coulomb potential V p for protons. Besides V^ p \p,/3) and V p the 
nucleon mean-held V n ^ p \p 1 (5) also includes a symmetric term for which we use a Skyrrne 
parameterization, i.e., 


= a(p/p „) + b(p/p 0 y +V?+ VZP(p.P). 


( 6 ) 


In the above, the parameters a, b and a are determined by the saturation properties and the 


compressibility K of symmetric nuclear matter |HJ . The symmetric term should also contain 
a momentum-dependent part. However, it is not essential for present study as we will show 
that the ratio of preequilibrium neutrons to protons from collisions of neutron-rich nuclei is 
essentially independent of the symmetric part of the nuclear EOS. The isospin-dependent 
BUU model was used to explain successfully several isospin-dependent phenomena in heavy- 
ion collisions at intermediate energies |35|. More recently, the isospin-dependence of collec¬ 
tive how and balance energy predicted in ref. [[5(J using this model was confirmed experi¬ 
mentally at NSCL/MSU by Pak et al. |37| , |38| . 


We have studied collisions of ll2 Sn + 112 Sn, 124 Sn + 124 Sn and 132 Sn + 132 Sn reactions 
at a beam energy of 40 MeV/nucleon. The hrst two reactions have been recently studied 
experimentally at NSCL/MSU by the MSU-Rochester-Washington-Wisconsin collaboration 
Preequilibrium particles were measured in these experiments and are now being an¬ 
alyzed m ■ The last system is included only for the purpose of discussions and comparisons. 
To identify free nucleons from those in clusters, we use in our calculations a phase-space coa¬ 
lescence method at 200 fm/c after the initial contact of the two nuclei, when the quadrupole 
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moment of the nucleon momentum distribution in the heavy residue is almost zero, indicat¬ 
ing the approach of thermal equilibrium. A nucleon is considered as free if it is not correlated 
with other nucleons within a spacial distance of A r = 3 fm and a momentum distance of 
A p = 300 MeV/c. We have checked that the results are not sensitive to these parameters if 
they are varied by less than 30% around the above values. 

We first study effects of the compressibility K of symmetric nuclear matter and the in¬ 
medium nucleon-nucleon cross section on the ratio R n / P (E]dn) by dropping both the Coulomb 
and symmetry potentials in the BUU model. In Fig. this ratio is shown as a function of 
nucleon kinetic energy for central (upper window) and peripheral (lower window) collisions 
of 132 Sn + 132 Sn at a beam energy of 40 Mev/nucleon. When varying the compressibility K 
from 210 MeV (open squares) to 380 MeV (filled circles), we find that although the yields 
of both neutrons and protons increase, their ratio remains almost the same for all impact 
parameters. This is simply because the effects of symmetric EOS on both neutrons and 
protons are identical. 

The experimental cross section for neutron-proton collisions is about three times that for 
neutron-neutron (proton-proton) collisions in the energy range studied here. Setting the two 
cross sections equal (fancy squares), we find that the yields and their ratios change by less 
than 10% even in peripheral collisions of 132 Sn+ 132 Sn. This result is also easy to understand 
since both colliding nucleons have the same probability to gain enough energy to become 
unbound (|0|. Thus, the in-medium, isospin-dependent nucleon-nucleon cross sections do 
not affect much the ratio R. n /p(Ekm) ■ It is important to point out that in the absence of 
Coulomb and symmetry potentials the ratios are almost independent of the nucleon kinetic 
energy and have a constant value of about 2.1 ±0.3 in both central and peripheral collisions 
of 132 Sn ± 132 Sn. 

Including the Coulomb and the asymmetric term of the EOS in eq. (|6|), one can then 
study the effects of the symmetry energy S(p) since the Coulomb effect is well-known. 
We expect that the symmetry potential will have the following effects on preequilibrium 
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nucleons. First, the symmetry potential tends to make more neutrons than protons 

unbound. One therefore expects that a stronger symmetry potential leads to a larger ratio 
of free neutrons to protons. Secondly, if both neutrons and protons are already free, the 
symmetry potential makes neutrons more energetic than protons. These effects are shown 
in Fig. m where we display the ratios calculated using the three forms of F(u) for central 
(left windows) and peripheral (right windows) collisions of 112 Sn + 112 Sn, 124 Sn + 124 Sn 
and 132 Sn + 132 Sn, respectively. The ratios change continuously from central to peripheral 
collisions, more detailed study on the impact parameter dependence of the ratio R n / P (E k[n ) 
will be published elsewhere |IT|. The increase of the ratios at lower kinetic energies in all 
cases is due to Coulomb repulsion which shifts protons from lower to higher kinetic energies. 
On the other hand, the different ratios calculated using different F(m)’s reflect clearly the 
effect mentioned above, i.e., with a stronger symmetry potential the ratio of preequilibrium 
neutrons to protons becomes larger for more neutron rich systems. 

It is interesting to note that effects due to different symmetry potentials are seen in 
different kinetic energy regions for central and peripheral collisions. In central collisions, 
effects of the symmetry potential are most prominent at higher kinetic energies. This is 
because most of the finally observed free neutrons and protons are already unbound in the 
early stage of the reaction as a result of violent nucleon-nucleon collisions. The symmetry 
potential thus mainly affects the nucleon energy spectra by shifting more neutrons to higher 
kinetic energies with respect to protons. In peripheral collisions, however, there are fewer 
nucleon-nucleon collisions; whether a nucleon can become unbound depends strongly on 
its potential energy. With a stronger symmetry potential more neutrons (protons) become 
unbound (bound) as a result of a stronger symmetry potential, but they generally have 
smaller kinetic energies. Therefore, in peripheral collisions effects of the symmetry potential 
show up chiefly at lower kinetic energies. For the more neutron-rich systems the effects of the 
symmetry potential are so strong that in central (peripheral) collisions different forms of F(u) 
can be clearly distinguished from the ratio of preequilibrium neutrons to protons at higher 
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(lower) kinetic energies. However, because of energy thresholds of detectors, it is difficult to 
measure low energy nucleons, especially neutrons. Furthermore, the low energy spectrum 
also has appreciable contribution from equilibrium emissions. Therefore, the measurement 
of the ratio i? n /p(^km) in neutron-rich, central heavy-ion collisions for nucleons with energies 
higher than about 20 MeV is more suitable for extracting the EOS of asymmetric nuclear 
matter. 

In conclusion, collisions of neutron-rich nuclei at intermediate energies reveal novel in¬ 
formation about the EOS of asymmetric nuclear matter that is of interest to both nuclcar- 
and astro-physics, such as the properties of radioactive nuclei, supernovae, and neutron 
stars. The ratio of the number of preequilibrium neutrons to that of protons is found to be 
sensitive to the asymmetric part, but not to the symmetric part of the nuclear EOS. It is 
also almost independent of the in-medium nucleon-nucleon cross sections. The experimental 
measurement of this ratio therefore provides a valuable means for determining the EOS of 
asymmetric nuclear matter. 
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FIGURES 



FIG. 1. The symmetry potential for neutrons and protons corresponding to the three forms of 
F(u) (see text). 
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FIG. 2. The ratio of preequilibrium neutrons to protons as a function of nucleon kinetic en¬ 
ergy for central (upper window) and peripheral (lower window) collisions calculated without the 
Coulomb and symmetry potentials. 
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FIG. 3. Same as in Fig. |] but calculated with both the Coulomb and symmetry potentials. 
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